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We study a nonequilibrium equation of states of an ideal quantum gas confined in the cavity under 
a moving piston with a small but finite velocity in the case that the cavity wall suddenly begins to 
move at time origin. Confining to the thermally-isolated process, quantum non-adiabatic (QNA) 
contribution to Poisson's adiabatic equations and to Bernoulli's formula which bridges the pressure 
and internal energy is elucidated. We carry out a statistical mean of the non-adiabatic (time- 
reversal-symmetric) force operator found in our preceding paper (K. Nakamura et al, Phys. Rev. 
E83, 041133 (2011)) in both the low-temperature quantum-mechanical and high temperature quasi- 
classical regimes. The QNA contribution, which is proportional to square of the piston's velocity 
and to inverse of the longitudinal size of the cavity, has a coefficient dependent on temperature, gas 
density and dimensionality of the cavity. The investigation is done for a unidirectionally-expanding 
3-d rectangular parallelepiped cavity as well as its 1-d version. Its relevance in a realistic nano-scale 
heat engine is discussed. 
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PACS numbers: 05.30.-d, 05.70.Ln, 51.30.+i. 
I. INTRODUCTION 

The equation of states plays an important role in ther- 
modynamics and statistical mechanics. Let's consider 
the Carnot's thermodynamic cycle proposed almost two 
centuries ago [l| . It is the most efficient cycle for convert- 
ing a given heat into work. In this cycle, the system is 
assumed to undergo a series of different thermodynamic 
states and performs work on its surroundings, thereby 
acting as a Carnot heat engine. However, such a perfect 
engine is only a theoretical limit and practical engines 
must incorporate the effect of non-zero velocity of the 
moving piston. 

In the Carnot cycle, the pressure (P) and volume (V) 
of an ideal classical gas (Boltzmann gas) confined in 
the cavity are assumed to obey the equilibrium equa- 
tion of states, i.e., Boyle- Charles' law (BCL) and a set 
of Poisson's adiabatic equations in the isothermal and 
thermally adiabatic processes, respectively. The Pois- 
son's adiabatic equations are derived from the first law 
of thermodynamics together with BCL. BCL itself is a 
special limit of the Bernoulli's formula (BF) bridging 
between the pressure (P) and internal energy (U) for 



quantum and classical gas in the cavity in d-dimensions. 
BF is available from the relation PV = — f2 with use of 
density of states in calculating the thermodynamic po- 
tential Q for both classical and quantum gas. To be 
specific, PV = \U, U, and 2U for d=3,2, and 1, re- 
spectively. The last case may be better rewritten as 
FL = 2U with use of the force (F) and the length (L) 
of the 1-d cavity. For a classical gas, U = ^NkT, NkT 
and TjNkT for d=3,2 and 1, respectively, with use of 
the number of particles N, Boltzmann constant k and 
temperature T. Then the Bernoulli's formula reduces to 
BCL, PV = NkT, irrespective of dimensionality For 
a quantum gas, the Bernoulli's formula works as well, 
where U = E (l + 0.0713 (mT/h 2 ) 2 (V/iV) 4/3 ) with 

E = (3/10) (6tt 2 ) 2/3 (h 2 /m) (N/V) 2 ^ 3 N for d=3 Fermi 
gas in the low-temperature and high-density regime (see 
Landau-Lifshitz j2|). In the thermally adiabatic pro- 
cess, a set of Poisson's adiabatic equations also works, 
which are given by PV r ( d+2 '/ <i =const., T(d ^ 2)/2 =const., 
and PT d / 2 =const., irrespective of classical and quantal 



systems [3 

In constructing Bernoulli's formula, the velocity of the 
wall of a gas container (cylinder, cavity, billiard, etc.) is 
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assumed to be negligibly small. To make the theory of 
heat engines more realistic, one must evaluate the effect 
of the non-zero velocity of the piston, i.e., the wall motion 
of the gas container. Since the kinetic theory of Boltz- 
mann gas tells that a moving piston does not play a role 
in the equation of states, we shall investigate the nona- 
diabatic dynamics in the quantum heat engine. While 
in recent years there appeared papers which treated the 
quantum engine, they were either concerned with a quan- 
tum analog of Carnot's engine [i|~0| or with a quantum 
analog of nonequilibrium work relation (i.e., fluctuation 
theorem) 0, 0| . And no work so far was engaged in nona- 
diabatic force and pressure due to a moving piston and in 
the statistical treatment of a noninteracting Fermi gas. 

In this paper, confining ourselves to the thermally- 
isolated process, we shall investigate the non-equilibrium 
equation of states for an ideal quantum gas (Fermi 
gas) confined into an expanding cavity in the case that 
the cavity wall suddenly begins to move at time ori- 
gin. Quantum non-adiabatic (QNA) contributions to the 
Bernoulli's formula and to Poisson's adiabatic equations 
due to the non-zero velocity of the moving piston is eluci- 
dated. In Sec. HIl with use of the nonadiabatic force op- 
erator in our preceding paper the adiabatic and non- 
adiabatic pressures are defined. In Sec lIIII expectation of 
nonequilibrium pressure is expressed in terms of density 
of states, which will enable the calculation of thermo- 
dynamic and statistical averages of nonadiabatic pres- 
sure. In Secs lIVI and [Vj the unidirectionally-expanding 
cavities in d=l and d=3 dimensions are studied and ex- 
plicit forms for QNA contributions to the pressure, in- 
ternal energy and equation of states will be given in the 
low-tempcrature and high density regime as well as in 
the high-temperature and low-density regime. In Sec lVIl 
physical implications of QNA contributions will given. 
Section IVIII is devoted to summary and discussions. In 
Appendices lAl and lEl using the Fermi-Dirac distribution, 
we summarize several formulas for thermodynamic aver- 
ages in the cases of 1-d and 3-d rectangular cavities. 

The isothermal process which requires a contact with 
the heat reservoir is outside the scope of the paper and 
will be investigated in due course. 



II. ADIABATIC AND NONADIABATIC 
PRESSURES 

Before embarking upon the adiabatic and nonadiabatic 
pressures, we shall briefly summarize the derivation of the 



adiabatic and nonadiabatic force operators in our preced- 
ing paper [10], but here in the context of the parallel- 
piped rectangular 3-d cavity. Let's consider a Fermi 
gas (non-interacting Fermi particles) confined in a cavity 
with a moving wall (i.e., piston). The wall receives the 
force from the Fermi gas in the cavity. Under the condi- 
tion that whole system consisting of Fermi particles and 
a moving wall keeps the energy conservation, the work 
done on the wall by the force is supplied by the excess 
energy due to the energy loss of Fermi particles showing 
the non-adiabatic transition. In this way one can con- 
ceive both the adiabatic and nonadiabatic forces. In the 
adiabatic limit, the adiabatic force due to the quantal gas 
on the cavity wall is proportional to the derivative of the 
confining energy with respect to the cavity size. What is 
a characteristic feature of the nonadiabatic force coming 
from the non-adiabatic transition? 

We choose a 3-d rectangular parallelepiped cavity with 
the size L% x Ly x L z 7 

one of whose walls is moving in 

x-direction (see Fig. [T]). 




FIG. 1: 3-d rectangular parallelepiped cavity confining the 
quantal gas, with the size L x (t) x L y x L z . One of its walls 
is moving in z-direction. P x (t) and F x (t) stand for the x- 
components of pressure and force. 

The original Hamiltonian for the cavity with a time- 
dependent longitudinal size L x {t) is given by 



Htotal = H + H± 



(1) 



with 



2m dx 2 1 

2m dy 2 dz 2 

The wavefunction is a product of the longitudinal and 
perpendicular parts: 



H 



(2) 



(3) 
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which satisfies the moving and static Dirichlet boundary 
conditions for -0 and ip±, respectively as 

i/)(x = 0,t) = i/>(x = L ai (t),t)=0, (4) 
^x(y = 0,z,t) = ipx(y = L y ,z,t) = 0, (5) 
i; ± (y,z = Q,t) = ip x (y,z = L z ,t) = Q. (6) 

Throughout the time evolution, the instantaneous (adi- 
abatic) eigenstate is characterized by a set of quantum 
numbers (n x ,n y ,n z ). The longitudinal perturbation in 
H commutes with the perpendicular part H± in the total 
Hamiltonian in Eq.([T|), and thereby the quantum num- 
bers n y , and n z are conserved against an expansion along 
x. Therefore, if a confined particle is initially in a mani- 
fold with the fixed n y and n z and the cavity expands only 
in x-direction, there occurs no mixing among manifolds 



with different 



and 



Consequently the dynamics 



of iptotai(x,y, z,t) is determined by the time-dependent 
Schrodinger equation for the longitudinal part ip{x,t) as 



,8 ,. . 
in—ip(x, t) 



Hip(x,t). 



(J) 



The expectation of the normal component of the force 
acting on the wall is obtained by 



F x 







dL x (t) 



{i>m) = -{i>\ 



dH 

dL x (t) 



M, (8) 



where, in obtaining the last expression, we used -gj^-\4 > ) = 
T~~nt\^) = ~TF~ HW) an d its Hermitian conjugate. 
Hence the force operator is defined by 



F x = - 



dH 

dL x {t\ 



(9) 



Since the original Hamiltonian H for the cavity with its 
time-dependent longitudinal size L x (t) does not formally 
include L x (t) explicitly, however, there is no way to define 
the force operator directly by using Eq.©. 

To overcome this difficulty, we shall make the time- 
dependent canonical transformation of H related to the 
scale transformation of both the coordinate x and am- 
plitude of the wave function ip. This transformation is 
defined bv[ll| 



H 



e -iV {H _ ih 9_ )e iU j 
px)hxL x {t) = i (x-S- 



(10) 



with U = -^(xp + px) hiL x (t) = i [x-§^ + ±) \nL x (t). 
This canonical transformation leads to the scaled coordi- 
nate x defined by e~' lU xe' lU = xL x (t) and the scaled wave 



function cj>(x,t) 



- -iU 



4>{x,t) = ^/L^^}{xL x ,t), The 



range of x is < x < 1, which is time- independent. 



Also the normalization factor of 4>(x,t) becomes L x - 
independent and satisfies the fixed Dirichlet boundary 
condition 0(0, t) = = 0. 

Finally the Schrodinger equation is transformed to 



iH d 4 = H4> 
ot 



with the new Hamiltonian 



H = 



2 



ImUi dx 2 



■^L x _ d ih L x 
L x dx 2 L x 



(11) 



(12) 



which is Hermitian. (f>(x, t) now satisfies the fixed Dirich- 
let boundary condition 0(0, t) = 0(1, t) = 0. 

Taking L x derivative of H, we can rigorously define 
the force operator in the transformed space, whose in- 
verse canonical transformation gives the force operator 
expressed in the original space as 



F 



P 

mL 



2L 2 



{xp + px) 



■--i 
m L 



i%k^(xd x + d x x), 



2L 2 



(13) 



where we suppressed the suffix x in both the force op- 
erator and the longitudinal length. The issue in Eq.(fl3|) 
is universal, irrespective of the kind of canonical trans- 
formations. In fact, one may choose another canonical 
transformation such as a combination of U in Eqs.(|10p 
and the gauge transformation (see [10(), which also guar- 
antees the wave function to satisfy the fixed Dirichlet 
boundary condition and the transformed Hamiltonian, 
say H', to be Hermitian. The derivative of H' w.r.t. L x 
defines the force operator F' , and the inverse of a combi- 
nation of the gauge and scale transformations results in 
the identical expression for F. 

In the final expression of Eq. flTB")) . the first and the sec- 
ond parts defines the adiabatic and nonadiabatic forces, 
respectively. The latter part, which gives an essential 
contribution when the system is not in the instantaneous 
eigenstates, is invariant under the time-reversal operation 
since both L and p change their signs. The expression 
in Eq.(|13[) is the force normal to the wall, and, when di- 
vided by an area of the wall, it gives the adiabatic and 
nonadiabatic pressures (P) acting on the moving wall of 
the 3-d rectangular parallelepiped cavity: 



P 



LyL z 



(14) 
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III. EXPECTATION OF NONEQUILIBRIUM 
PRESSURE IN TERMS OF DENSITY OF STATES 

Let' s consider the system to be thermally isolated and 
the wall of the cavity to begin to move at t = suddenly 
(see Fig. The Fermi gas in the cavity is assumed 

to satisfy the equilibrium Fermi-Dirac distribution until 
t = 0. The expectation of the force operator is evaluated 
in terms of the density operator p: 



F = Tr(pF). 



The density operator p for a thermally-isolated 
nonequilibrium state of the Fermi gas obeys von Neu- 
mann equation 

dp 



ifj- 



M 



[H,P] 



(16) 



where the original Hamiltonian H and coordinate x are 
used. 



L x (t)k 



L x (0) 



L x (0)+L x t 







FIG. 2: Time dependence of L x (t). The wall is assumed to 
begin to move at time origin. 

Using adiabatic bases (instantaneous eigenstates) 
{|n)}, the matrix elements of p satisfies 



— Jfr(E n — E m )p nm — — | ^2 llnPlm + ^ 



With lmn = (-l)m+n+l 2mn (J _ ^ 



Then F becomes 



E ^ ^ Pnn\En 



where 



h 2 (mr) 2 ihL 

^ ran — ®n\n ~T T 'Jmn • 



(17) 



(18) 



(19) 



To make the problem tractable, we assume L < dj?, 
that is, the wall velocity L is much less than the Fermi ve- 
locity vf , which guarantees a confined particle to collide 
with the cavity wall many times during the wall displace- 
ment of O(L). The above unequality is scaled by L and 
is written as 



L 1 
- « — 

L tf 



(20) 



(15) where Tp(= — ) is a characteristic time for a particle 



to travel through the cavity. With use of the smallness 
parameter ^, we substitute the expansion 



p = f(H) + - 9l 



L 



<I2 



(21) 



into the von Neumann equation. Then, for orders of (f;) , 
Y and (£) 2 , we have f(H) = 0, g lnm = E \ h E ™ g lnm - 

{imnfm + lnm.fn), and g 2 nn = - J2 ll 'tnidUn + 9lni), 

respectively. 

The condition Eq. (|2T)]) guarantees Lt <C L in a wide 
time range. Then a set of the above equations can be 
solved as, 



1 



s /3(E„-/i) _|_ I 



8ntr 



fn$n 



(22) 



Q lmn, — 



92, 



E n — E n 



(l - e^^tj (/„ - f m ), (23) 



= -2x;-&(/»-/*) 

E n — E( 
x 1 — cos 1 



E n — Ei 



(24) 



Equation (|22[) denotes the initial Fermi-Dirac distribu- 
tion with inverse temperature P = -gj:- 
Expectation value F is given by 



F = JFi + F 2 + F3 



(25) 



where 



5 



Afo2 ^n-^m fn frn 

£3 (E n - E m f E n - E v 



L 3 



ra^n 



- Tp p 



(27) 



ft 2 / L 



f 3 - - (jj E^nn^n 



, , > / -En E m 

x 2 sin 1 

2/i 



-327r 2 ft 2 ^ £ E n E m ^—^-6iE„ E,„ i. 



(28) 



It should be noted that both and F3 are quadratic in 
L. The absence of L-linear terms is caused by a subtle 
cancellation of the linear cross-coupling terms among the 
matrix elements of the force operator and those of the 
density matrix both expressed as a series expansion w.r.t. 
L. 

In obtaining the final expression for Fx, F2 and F3, we 
used E n = 1T 2 mL' 1 anc ^ * ne as y m Pt°tic form, 



AE 



n5{AE), 



(29) 



which is valid in the time domain much larger than 
the minimum resolution of time (t S> xt?)' Thanks to 
Eq. (f2U)) . the explicit time dependence of F 2 and F3 are 
suppressed. 

The discrete summations can now be reduced to con- 
tinuum integrations with use of 1-d density of states as 



E = E 



D 1 (E)dE. 



(30) 



Noting E 



£>i (E) is given by 



Using the above facts, we shall write the final results for 
A, F 2 and F 3 : 

2 



Ft 



EDi(E)f(E), 



(32) 



F 2 = -8ir 2 h 2 — / dE / dE'EE' -f- 
L d Jo Jo dE 

x D l (E)D 1 (E')8{E - E'), 



E=E' 



(33) 



-32tt 2 /i 2 



I? 



dE 



dE'EE' 



dE 



E=E' 



x D 1 (E)D 1 (E')8(E - E') 



(34) 



The purpose of the present paper is to generalize 
Bernoulli' s formula bridging pressure and internal en- 
ergy to the case of the expanding cavity. Therefore, one 
should also provide general formulas for the internal en- 
ergy in the case of a moving piston. With use of the 
expansion in Eq. (|21|) and the matrix elements 



(H)r 



h 2 n 2 n 2 



E n 5 n 



2mL 2 n 
we have the internal energy U, 

U = Tr(pH) = U 1 + U 2 + U 3 . 

Here 



(35) 



(36) 



OO />OG 

Ui = Tr(fH) = E E n f n = / ED l {E)f{E)dE (37) 



and 



L 



L 1 E n — E v 

n>m 



,IE I dE'EE' ^- 
dE 



x D 1 (E)D 1 (E')S(E- E'). 
Noting the absence of the diagonal elements of g\ , 

U 2 = jTr( 9l H) = 0, 
namely, always vanishing. 



E=E> 

(38) 



(39) 



IV. CASE OF EXPANDING 1-D CAVITY 

Firstly, concentrating on the expanding 1-d cavity, we 
shall evaluate the final expressions in the previous section 
in two limiting cases, i.e., in the low-temperature and 
high-density region for a degenerate quantum gas and in 
the high-temperture and low-density region for a quasi- 
classical gas. 
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A. Low-temperature and high-density region 



and 



Having recourse to formulas in Eqs. (|Al| and (|A2[) . the 
expectation of force terms in Eqs. ([32j) - (j34|) becomes 



Fi = 



2Jm 



^/2nh Jo 
2v/2m 



3nh 



-/' 



3/2 



Ef(E)dE 

^2 



7tt 4 



1 + '—{kTf^ 2 + ^-(fcT)V^ 4 + 



640 



we have 

V2mL 



3nh 



-J" 



3/2 



(40) 



Us = 2mL ix, 



-{kTfix- 



(47) 



2mL 2 /x. 
(48) 



F 2 = —4m- 



lr 



df 

dEE— 
dE 



Fa 



-8m- 



Zr 



dEE 



df_ 

' dE 



L 2 
4m— /i, 



5m— /x. 
L 



(41) 



(42) 



Using the expansion for /x in Eq. (IA4p , we have 



17 = 



ir 2 h 2 
6m 



AT 1 + 



ir 2 H 2 L 2 I ^ 



where \x is the chemical potential. Thereby, 




4- 12m— 
L 



2 v / 2m 
37rft ' 



,3/2 



1 + — (fcT)V 



(49) 



The first term corresponds to the 1-d version of the ex- 
(43) isting result (Landau-Lifshitz[2j), and the second one is 
a nonequilibrium correction. Combining Eqs. (|44[) and 
Noting the low-temperature and high-density expansion (|49|) . we have 
of [x with use of particle number N in Eq. (|A4[) , we have 

- 2U 

ATr 2 h 2 L 2 ' 



F = 




1 fmkTY 



1 



1 

3k 2 



mkT 



r, 2 



Equation (|4"4")) can be rewritten as 



FL 



ir 2 H 2 
3m 



-N 



l + -(^ 



^44^which generalize the Bernoulli's formula in 1-dimension. 
The right-hand side gives a qualitatively new correction 
term due to a moving piston. This equation stands for 
the nonequilibrium equation of states for a quantal gas 
confined in the expanding cavity with the finite velocity 
(L) of a piston. 



6Tr 2 h 2 L 2 N 2 



1 + 



kTY 
12 \~) 



(45) 



Noting — =const. in the unperturbed adiabatic state, 
Eq.([4"5j) is nothing but a generalization of Poisson's adi- 
abatic equation (PAE) in 1 dimension in the case that a 
piston has a small but non-zero velocity. A qualitatively 
new correction term on the right-hand side is quadratic 
in both velocity of the piston and particle number. 

The internal energy for the expanding 1-d cavity is 
calculated in a similar way: Noting 



B. High-temperature and low-density region 

In this subsection we shall investigate the opposite 
limit, i.e., the high-temperature and low-density quasi- 
classical regime. 

Here we shall have recourse to a high-temperature 
expansion of Fermi-Dirac distribution expansion with a 
negative value fx, 



Ui = 



/2mL 
3irh 



3/2 



1 



-(kTfii- 



(46) 



f(E) 



1 



j = 5Z(-l) n_1 e"" 



(51) 



(50) 
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Substituting Eq. (|5"Tj) into middle terms in each of 
Eqs. p0|) - (|42p . one can evaluate the force: 



V2nh^ Jo 



2Trh 2 



(kT)f e P» 1 



/2 e -n0(E-ri dE 



2V2 



0{e 



(52) 



-nP(E-fj,) 



f 2 00 /-do 

F2+F3 = -12m— V(-l) n n/3 / dESe 
L n=l ' /o 

= Um^kTe^ (l - ^ + 0(e 2 ^)") . (53) 



Therefore 

-F = Pi + F\ 2 + -F 3 



(^)< tr » 3/v "( 1 -|J + ^") 

L 2 



12m- 



2V2 

-fcTe^ C 1 - -e^ + 0(e 2/3 ") 
i/ V 2 

Equation (p3|) can be rewritten as 



(54) 



{kTf/ 2 



e^(l 



2^ 2 y 

L 2 



12 TO _(fcT)- 1 / 2 e^(l 
L 



)• (55) 



Since /3/x =const. in the unperturbed adiabatic state 
(see Landau-Lifshitz [2]), Eq. (j55|) is a generalization the 
Poisson equation in 1 dimension expressed in terms of 
pressure and temperature. 

Using in Eq.(f5~3| a high-temperature and low-density 
expansion of e' 3 ' 1 in Eq. (|B3[) . we have 




[2 N 

Yl^Phih VmkT 

L L \ 



\p2~ ) V mkT J 
(56) 



Similarly, the internal energy is now given by 



N ( 3^ f% 

= —kT 1 - ^ 

2 \ 2 v mfcT 



2V2^hL 2 — VmkT 1 - 1 _ , N ., 




(57) 



Therefore, a generalized Bernoulli's formula in the 
quasi-classical region is given by 



,N 



1 



FL-2U = 8V2^HL 2 — VmkT 1 - ( 1 + -= ] y^F ; L 

i \ V V2/ V / mjfeT / 

(58) 

The right-hand side is a nonequilibrium contribution due 
to the finite velocity of a piston. We find that a deviation 
from Bernoulli's formula appears only when the quantum 
effect will be incorporated. In fact, in the limit h 0, 
we see U = ^kT and Eq. ([55|) becomes the 1-d version 
of Boyle-Charles' law, FL = NkT, which includes no 
contribution due to kinematics of the piston. 



V. CASE OF 3-D RECTANGULAR 
PARALLELEPIPED CAVITY SHOWING A 
UNIDIRECTIONAL EXPANSION 

The realistic heat engine is composed of a 3-d cavity 
with a piston moving in a fixed (x) direction. The force 
F in the previous sections is taken as x component of the 
force vector for the case of 1-d motion of the piston in the 
3-d rectangular parallelepiped cavity with size L x ^ Ly x 
L z under the fixed perpendicular (or transverse) modes 
(n y ,n z ). 

We shall denote F x as the x-component of the force 
vector averaged over both longitudinal and perpendicular 
modes. Noting Eq.(fT4|). the expectation of pressure on 
the wall of a piston is given by 



LyL z 



(59) 



where L y L z is an area of the wall. 

F x can be evaluated in a similar way as F, but Fermi- 
Dirac distribution should include a contribution of the 
energy due to the perpendicular modes. Namely, the 
eigen-energy of a particle is now 

E(n x , n y , n z ) = E\\ (n x ) + Ej_(n y ,n z ) (60) 



with 



E\\{n x ) 



fc2 / \ 2 



2m V L 



Ej_(n y ,n z ) 




(61) 



(62) 
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and Fermi-Dirac distribution is given by 

1 



f(E) = 



(63) 



The statistical average is the one over the longitudinal 
mode (n x ), followed by another one over the perpendic- 
ular modes (n y ,n z ). The expectation value F x is given 
by 



F = F i 4- F 

± X M Xl ~ 1 X 



F, 



x3 • 



(64) 



In the low-temperature and high-density regime, we 
have the following results: 

r-E 



This is a 3-d version of Poisson's adiabatic equation 
which now incorporates the non-adiabtic contribution. 
As in the case of the 1-d cavity the qualitatively new 
contribution is proportional to the square of the wall ve- 
locity and to inverse of the longitudinal size of the cavity, 
but the coefficient shows a different dependence on par- 
ticle number. 

The internal energy for the 3-d rectangular cavity with 
a moving piston is straightforward: 
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where we employed the 2-d density of states 
D 2 (E) 



E=E' n +E ± 

x D 1 (E ll )D 1 (E\ l )D 2 (E - E^SiE'w — £ , y).(72) 

It should be noted that, in the calculation of £/ 3 ~ d , the 
(66)bulk energy E(= E\\ + E±) is averaged which is a 3-d 
generalization of the 1-d energy. The final result for the 
internal energy is 



2L y L z m 



(67) 



7T ft 2 

together with D\(E) in Eq. (|3"T1) . F x ^ can be obtained in 
a similar way, but 8tt 2 H 2 j Q dE'^ in the integral of F x2 is 
to be replaced by 327r 2 ft 2 J*^" dE'^ , which eventually leads 
to F X 3 = 2F x2 . 

Then the pressure on the wall is 
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(With use of the expansion for fi in Eq. (|A6|) . the first 
term on r.h.s. of Eq. (j?3")) proves to agree with the result 
by Landau-Lifshitz reproduced in Introduction. The mi- 
nor discrepancy of a numerical pref actor of 0(1) is due 
to our choice of anisotropic density of states in Eq. (|7Tj) 
for the 3-d rectangular parallelepiped cavity.) 

The Bernoulli's formula in the present case becomes: 
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With use of the low-temperature expansion of [i in 
Ea.(|A6|l. Eq. (|551) can be written as 
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with V = L x L y L z . With use of the low-temperature 
expansion of /i in Eq. (|A6j) . Eq.(|74]) can be rewritten as 
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with cq = ^ „ and c\ 



The right- 



3 ^ 3(3V2) 4 / 3 ^ 2 / 3 ' 

hand side is a nonadiabatic contribution to the equilib- 
rium equation of states in 3-dimensions due to a moving 
piston. In the quantum adiabatic limit L x = 0, the above 
equation reduces to the standard Bernoulli's formula for 
the 3-d quantum gas. 

We shall proceed to the high-temperature and low- 
density regime. The values of F x i, F x2 , F x3 , U 3 ~ d and 
U 3 ~ d for the 3-d cavity are evaluated, by substituting 
into the second expressions in each of Eqs. (|65l) . AMI) . 
(fTTj) and (j?2"|) the expansion of Fermi-Dirac distribution 
in Eq. (j5l"j) . The results are: 
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Then the pressure defined by 
F x i + F x 
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Similarly, we see 
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leading to the internal energy, U 3 d = 
Bernoulli's formula is now given by 
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Noting the high-temperature and low-density expansion 
of in Eq. (|B5[) . we see 
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4\/2 V" \mkT / 
We can confirm that the nonadiabatic contribution (NC) 
appears as a quantum effect and plays a role with decreas- 
ing the system's size (L x ). In other words, NC vanishes 
in the classical limit (ft — >■ 0), which is consistent with the 
kinetic theory of Boltzmann gas which incorporates the 
effect of moving piston. The essential results obtained in 
this Section is summarized in Table 1. 



Table 1. Non-adiabatic contributions to equation of states in thermally-isolated process in 3 dimensions. 
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high-temperature quasi-classical region 
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VI. PHYSICAL IMPLICATIONS OF QUANTUM states in the cases of 3-d rectangular parallelepiped cav- 
NONADIABATIC CONTRIBUTIONS ity and its 1-d version, separately. To physically interpret 

the obtained results, however, it is more convenient to see 

So far we have obtained the completely analytical non- 
adiabatic contribution to the non-equilibrium equation of 
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the non-equilibrium equations of states for the general d- 
dimensional hyper-rectangular cavity which has the vol- 
ume V = L d ^ 1 L x and the moving wall (surface) with 
area S — U 1 ^ 1 . Such general derivation is also possible 
by using the density of states in d dimensions. Af- 
ter tedious and lengthy calculation (to be published else- 
where), the Bernoulli's formulas for the d-dimensional 
cavity are given by 



PV 



d L x 
Lj 



(84) 



and 



4) QNA contributions play an essential role in 

Bernoulli's formula rather than in Poisson's equation. In 

f 2 

fact, the coefficients prior to -£- are increased in Eq. 
(l84l and decreased in Eq. (l86l) as particle density 
is increased. Similarly, the coefficients is increased in 
Eq. ([55]) and decreased in Eq. (|87)l as temperature is 
increased. 

The above 4 issues constitute a punchline of the present 
paper. The Poisson's adiabatic equation and Bernoulli's 
formula, both of which are the basic laws of thermody- 
namics, are now generalized so as to include the QNA 
contributions that have never been reported so far. 



PV-~U~h(mkT)i(^)V^- ) (85) 

respectively, for the low-temperature high-density and 
high-temperature low-density regions. In a similar way, 
the corresponding Poisson's adiabatic equations are 
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and 
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respectively, for the low-temperature high-density and 
high-temperature low-density regions. The apparently- 
extra dimensionality of energy ([ML 2 T -2 ]) on the right- 
hand sides in all of the four equations above is traced 
back to our simplified replacement in Eq. (f2T)f and there- 
fore can be suppressed. Equations (|84 p -(|87 p recover all 
the results for d = 1 and d = 3 cavities in the previous 
Sections. We find the important features: 

1) Quantum non-adiabatic (QNA) contributions are 
quadratic in the wall velocity and therefore time-reversal 
symmetric, in marked contrast to the conventional be- 
lief 12j that the nonadiabatic force is linear in the wall 



velocity and breaks the time-reversal symmetry; 

2) QNA contributions are positive, which means that 
the moving wall gives rise to the apparently repulsive 
interaction among non-interacting Fermi particles, irre- 
spective of the direction of the wall motion, namely for 
both expansion and contraction of the cavity; 

3) QNA contributions are inversely proportional to the 
longitudinal size of the cavity and become more and more 
important when the cavity size is decreased. In partic- 
ular, they will play a nontrivial role in nano-scale heat 
engines based on quantum dots; 



VII. SUMMARY AND DISCUSSIONS 

Confining ourselves to the thermally-isolated process, 
we study a nonequilibrium equation of states of an ideal 
quantum gas confined to the cavity under a moving pis- 
ton with a small but finite velocity. The cavity wall is 
assumed to begin to move suddenly at time origin. Quan- 
tum non-adiabatic (QNA) contribution to Bernoulli's for- 
mula which bridges the pressure and internal energy is 
elucidated. Statistical means of the non-adiabatic (time- 
reversal symmetric) force and pressure operator [lo| 
are carried out in both the low-temperature quantum- 
mechanical and high temperature quasi-classical regimes. 
QNA contributions are quadratic in the piston's velocity 
and therefore time-reversal symmetric, in marked con- 



trast to the conventional belief 12], and they are posi- 
tive, which means that the moving piston gives rise to the 
apparently repulsive interaction among non-interacting 
Fermi particles, for both expansion and contraction of 
the cavity. QNA contributions are inversely proportional 
to the longitudinal size of the cavity, and thereby play a 
nontrivial role in nano-scale heat engines based on quan- 
tum dots. The investigation is done for an expanding 
3-d rectangular parallelepiped cavity as well as its 1-d 
version. The nonequilibrium contributions to Poisson's 
adiabatic equation are also elucidated. 

In the context of a classical gas, Curzon and 
Ahlborn 



13] and others 



14 



16j investigated a finite-time 



Carnot heat engine and obtained an interesting efficiency. 
However, they neither considered a quantum gas nor 
showed a nonequilibrium equation of states due to a mov- 
ing piston. Therefore one of the directions to extend our 
work may be to proceed to the same analyses as given 
here of the isothermal process which requires a contact of 
nano-scale engine with a heat reservoir. Another direc- 
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tion may be the fast-forwarding of the adiabatic expan- 
sion of a cavity [171 Il8l ] in the framework of von Neumann 
equation, to see an accelerated quantum Carnot heat en- 
gine. These subjects will be investigated in due course. 

Acknowledgments. One of the authors (K. N.) is grate- 
ful to Adolfo del Campo, Takaaki Monnai, Takahiro 
Sagawa and Ayumu Sugita for useful comments. 



which leads to the low-temperature expansion of [i as 
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Appendix A: Thermodynamic averages in 
low-temperature region at T <C To (degenerate 
temperature) 

With use of the Fermi-Dirac distribution f(E) = 
e <3(-E-M)+i > we summarize several formulas for thermody- 
namic averages (see Landau-LifshitzQ) in the cases of 
1-d and 3-d rectangular cavities. 

In low-temperature region at T <C T (degenerate tem- 
perature), we see 
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where g(Eo) = ip(Eo) = is assumed. 

Choosing 1-d density of states Di(E) as g(E), we have 
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N = J D 1 (E)f(E)dE 
from which the chemical potential is obtained as 
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This expansion is justified in the low-temperature and 
high-density regime. 

In the case of the 3-d rectangular cavity, 
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Appendix B: Thermodynamic averages at 
high-temperature region at T> To 

In the case of high-temperature region at T ^ To, we 
shall use a high-temperature expansion of Fermi-Dirac 
distribution with a negative value fi as given in Eg. (1511) . 
Then we see 

/>oo °° />oo 

/ g(E)f(E)dE^J2(- 1 ^ 1 g(E)e- n ^ E -^dE. 

J Eq n— 1 J Eo 

(Bl) 

Choosing 1-d density of states D\(E) as g(E), we have 
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from which fj, is determined by 
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This expansion is justified in the high-temperature and 
low-density regime. 

In the case of 3-d cavity, the particle number is 
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f^y\ e to (i _ , (B4) 



= V2^- 3 / 2 



and chemical potential is expanded as 
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